A unifying approach to left handed material design 
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In this letter we show that equivalent circuits offer a qualitative and even quantitative simple 
explanation for the behavior of various types of left-handed (or negative index) meta-materials. This 
allows us to optimize design features and parameters, while avoiding trial and error simulations or 
fabrications. In particular we apply this unifying circuit approach in accounting for the features 
and in optimizing the structure employing parallel metallic bars on the two sides of a dielectric film. 
Pacs: 42.70.Qs, 41.20.Jb, 42.25.Bs, 73.20.Mf 



Left-handed materials exhibit a negative permeabil- 
ity, /i, and permittivity, e, over a common frequency 
range Negative permeability is the result of a strong 
resonance response to an external magnetic field; neg- 
ative permittivity can appear either by a plasmonic or 
a resonance response (or both) to an external electric 
field. Negative [i and negative e lead to negative index 
of refraction, n, and to a left-handed triad of k, E, H; 
hence, the names negative index materials (NIMs) or 
Left-handed Materials (LHMs). Pendry^, Q suggested 
a double metallic split-ring resonator (SRR) design for 
negative /i and a parallel metallic wire periodic structure 
for an adjustable plasmonic response. Several variation 
of the initial design have been studied; among them a 
single ring resonator with several cuts has been proved 
capable of reaching negative /i at higher frequency Q ; in 
Fig. []Ja) a two cut single ring is shown schematically. 
This, by a continuous transformation, can be reduced to 
a pair of carefully aligned metal bars separated by a di- 
electric spacer of thickness t s HQ; in FigsCUb) and [He) 
the view in the (E, k) and (E, H) planes of this structure 
is shown together with the directions of k,E,H of the 
incoming EM field. 

The design shown in FigsQJb,c), besides its simplicity, 
has distinct advantages over conventional SRRs. The in- 
cident electromagnetic wave is normal to the structure 
as shown in Fig. HJb) , which enable us to build NIMs by 
only one layer of sample and achieve relatively strong re- 
sponse. Conventional SRRs, although they exhibit mag- 
netic resonance which may produce negative n, they fail 
to give negative e at the same frequency range and, hence, 
they are incapable by themselves to produce NIMs. An 
extra continuous wire is needed to obtain negative e via 
plasmonic response @,0|- In contrast, the pair of parallel 
metallic plates is expected to exhibit not only a magnetic 
resonance [Fig. Hfc), antisymmetric mode], but to show 
an electric resonance as well [symmetric mode] properly 
located in frequency by adjusting the length, I, of the 
pair. 

The simulations were done with the CST Microwave 




FIG. 1: (Color online) The two cut single metallic SRR (a) 
can be transformed to a pair of parallel metallic bars sepa- 
rated by a dielectric (b, view in (E, k) plane; c, view in (E, H) 
plane). By adding continuous wires, design d (view in (E, H) 
plane) results, which can be modified to a fully connected 
one on both sides of the thin dielectric board (e) . The dashed 
square defines the unit cells with dimension a x (parallel to 
H), a y (parallel to E) and a z (parallel to k). 



Studio (Computer Simulation Technology GmbH, Darm- 
stadt, Germany) using the lossy metal model for copper 
with a conductivity a = 5.8 x 10 7 for a single unit cell 
with periodic boundary in the (E,H) plane, field distri- 
bution and scattering amplitudes have been calculated. 
The e, \x in Fig. [6] have been obtained by a retrieval pro- 
cedure At the magnetic resonance the two plates 
sustain anti-parallel currents producing a magnetic field 
B confined mainly in the space between the plates and 
directed opposite to that shown in Fig. [He); the electric 
field, because of the opposite charges accumulated at the 
ends of the two plates, is expected to be confined within 
the space between the plates and near the end points. 
Indeed, detailed simulations, shown in Fig. HJc), confirm 
this picture. 

At the electric resonance the currents at the two bars 
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Notice that at telecommunication or optical frequen- 
cies, where the linear dimension are in the tens or hun- 
dreds of nm, the kinetic energy of the drifting electrons 
makes a contribution comparable or larger than the mag- 
netic energy. Hence, another additional inductance must 
be added to the right hand side of EqQ] Q • 

Each of the capacitance C m must be given by a formula 
of the type 



FIG. 2: (Color online) At the magnetic resonance the currents 
(a, in (E,k) plane, view in H direction), the magnetic field 
(b, in (H, k) plane) and the electric field (c, in (E, k) plane) 
are shown. Sizes of the cones show the intensity of magnetic 
field H (b) and electric field E (c) in logarithm scale. 



are parallel (symmetric mode); the magnetic field lines 
go around both bars, while the electric field is mostly 
confined in the space between the nearest neighbor edges 
of the two pairs of bars belonging to consecutive unit 
cells. 
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FIG. 3: (Color online) Current distribution of the two paral- 
lel metallic bar design (a) (side view, the parallel plates are 
behind one another) can be accounted for by the equivalent 
circuit (b), which, since points 1 and 2 are equivalent because 
of the periodicity, reduces to circuit (c) and (d) for the mag- 
netic (c) and electric (d) resonance respectively. 

The field and current configurations for both the an- 
tisymmetric and the symmetric mode can be accounted 
for by equivalent R, C, L circuits as shown in Fig. [3] (in 
which for simplicity the resistor elements have been omit- 
ted). Near the magnetic resonance frequency where the 
current configuration is as in Fig. [3jc), the magnetic field 
is between the two plates and it is, to a good approxima- 
tion, uniform (Fig. [2Kb)) . Hence the total inductance L, 
as calculated by the magnetic field energy, is 



L — 



t s , 
W 



(1) 



where I is the length of the wire, t s is the thickness of the 
dielectric spacer and w is the width of the wire. 
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where by inspection of Fig. [2fc), I — c\l with the nu- 
merical factor ci in the range 0.2 < c\ < 0.3. The capac- 
itance C e can be approximated by that of two parallel 
wires of radius t m and length w at a distance b apart 



C, = 



ln(6/*r 



(3) 



where t m is the thickness of each metallic bar and b is the 
separation of neighboring pairs Fig. [2Ka,b) (b = a y — I). 
The magnetic resonance frequency, uj m , is obtained by 
equating the impedance Z (of L m and C e in parallel) 
with minus the impendence —i/C m uj of the capacitance 

Cm- 

Since Z = iL m uj /(l — L m C e uj 2 ) we obtain 
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(4) 



The last relation follows because, for the values we have 
used (I = 7mm, w = 1mm, t s = 0.254mm, t m = 10/im 
and b = 0.3mm), C e ~ 0.1C m . Combining the EqQ] and 
Eql2]we find that 
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where e r = 2.53 is the reduced dielectric constant of 
the dielectric, e r = e/eo- I n Fig|4] we compare our re- 
sult of EqlH which shows that f m is a linear function of 
only l/l, with detailed simulations results. Fig. 2] shows 
the dependence of the magnetic resonance frequency as 
obtained from the retrieved resonant effective [i on the 
inverse length of the parallel metallic bars (Fig. [TJa), 
for different widths w (separation between parallel bars 
t s = 0.254mm fixed) and two different separations t s 
(width of bars w = 1mm fixed). Complete quantitative 
agreement is obtained if c\ = 0.22. Notice the indepen- 
dence of the simulation results on the width w and the 
dielectric thickness t s . It is worthwhile to point out that 
the result f m ~ l/l is robust over a wide range of pa- 
rameters even if Eqs. Q] and [2] are not valid. To see this 
point, consider the extreme case of a pair of thin wires 
(as opposed to a pair of bars) of length I, cross-section 
radius r, at a distance d apart such that r <C d <C I. For 
such a system L = (/z/47r)[l + 41n(d/r)] ~ ((iI/tt) ln(d/r) 
and C ~ eirl / \n(d / r) . Thus again f m ~ l/^fejll. 
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nance is much stronger and, hence, much wider we have 
to bring the magnetic resonance within the negative re- 
gion of e, i.e. we must have f e lower than f m as in Fig. 
IHIb), rather that the other way around, i.e. we must have 



L„ 



< 1 



(6) 



0.15 0.20 
1/1 (mm" 1 ) 

FIG. 4: (Color online) Linear dependence of the magnetic 
resonance frequency, f m , as obtained by simulation, on the 
inverse length I; this result as well as its independence on w 
and t s is in agreement with the simple formula {5). (t s — 
0.254mm for triangular, cross, circle; w = 1mm for diamond; 
and for all cases, b = 0.5 ~ 5.5mm, a x — 20mm). 



For frequencies near the electric resonance, because of 
mirror symmetry in Fig. [3fd), there is no current passing 
through the capacitances C m . As a results the electric 
resonance frequency f e is given by f e = \/(2ir^fC e L e ), 
where L e is expected to be of the form (fi/ir)g(w/l) where 
g{x) is a function which for x — > behaves as — ln(x). 

We point out that f e is a rather sensitive function of 
the small distance b, because C e depends on b, while 
L e is practically independent on b. Indeed the ratios 
f e (2b)/fe(b) and f e (3b)/f e (b) for b = 0.1mm according 
to the equation of C e and equation of f e are respectively 
1.14 and 1.215 in good agreement with the simulation 
results in Fig. [5] (1.13 and 1.21 respectively); the depen- 



This can be achieved by increasing C e either by decreas- 
ing b or by increasing at the ends of each bar the width 
w choosing a double T shape for each bar [9(. 

Still another possibility to make the negative e region 
wider (and more negative) is to add continuous metallic 
wires as in Fig. [IJd) which produce a plasmonic response 
[5J. By adjusting the width of these wires their effective 
plasma frequency f p can be made larger than the fre- 
quency, /i, at which the continuous curve in Fig. Etb) 
crosses the axis (/i ~ \QGHz). 
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and f e on a y /l = 1 + b/l is shown in points a (a y = 7.3mm 



FIG. 6: (Color online) Retrieved e c g (solid lines) and /u e ff 
(dotted lines) for two cut wires, (a) and (b) correspond to 
a x = 20mm) and b (a v = 7.1mm, 
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FIG. 5: (Color online) Magnetic resonant frequency f m cross 
over with electrical resonant frequency f e as a y /l — 1 + b/l 
varies between 7.1mm and 7.3mm; a x = 20mm. 

Fig. [5] in combination with Fig. [6] suggest the opti- 
mum design parameters for making the two bar scheme 
to produce negative index n: One has to avoid the cross- 
ing region where essentially to a considerably degree the 
two resonances cancel each other. Since the electric reso- 



a x = 20mm) in Fig. [S] Notice that both the response are 
Lorentz like. 

Finally the width of the bars, w, can increase until the 
bars join the "infinite" wires producing thus a continuous 
connected network which can be constructed by opening 
periodically placed rectangular holes on uniform metallic 
films covering both sides of a dielectric sheet [13, HH, E3] • 

In this letter we have shown that L, C equivalent cir- 
cuits can account for the EM properties of various neg- 
ative index artificial meta-materials (NIMs), even at a 
quantitative level; furthermore, this simple unifying cir- 
cuit approach offers a clear guidance in adjusting the 
design and optimizing the parameters for existing and 
possibly, future NIMs. 
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